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ABSTRACT 







The concept,, instrumentations and analytical techniques for the 
investigation of attenuation of wind generated water surface waves in 
a small scale real environment by separation of generation and propaga- 
tion by time in the dying wind situation were developed, and an empirical 
attenuation function applicable to a wide range of frequencies proposed. 
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I. INTRODUCTION 



The concept of wind-driven ocean waves as spectra of harmonic 
components which follow the laws of linear hydrodynamics was introduced 
in 1948 by Barber and Ursell (1), There has since been considerable 
progress toward the mathematical and physical description of ocean wave 
processes, with the ultimate objective being the prediction with reason- 
able certainty of future wave conditions at any geographical point on 
the basis of worldwide synoptic observations of both wave and wind fields. 
A summary of the major contributions to knowledge of ocean waves was 
compiled in 1967 by Cartwright (2). 

While theoretical models of wave generation, dispersion and geo- 
metrical spreading have been fairly well confirmed empirically and 
are used as the current basis of wave forecasting, the attenuation of 
waves during propagation is less well understood; and an attenuation 
function is consequently not applied in many wave —forecasting techni- 
ques . 

Some classical theories of wave attenuation, namely air resistance 
and turbulent friction, were critically examined by Darbyshire (4) and 
found wanting as singular explanations. The theory of non-linear scat- 
tering proposed by Phillips (6) appears promising and has been reason- 
ably supported by Snodgrass, et. al. (11) . Conceivably, the answer 
may lie in a combination of these processes. 

As has been the case in studies of other environmental processes, 
it appears that the phenomenon under consideration must be completely 
described empirically before the question of theory can be resolved. 

So far, environmental studies of attenuation have been confined to low 
frequency ocean swell, where the range of measurable frequencies is 
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included in only one of several logarithmic decades of possible surface 
wave frequencies. The classical theories and the generalized model 
developed later in this paper call for the attenuation function to be 
exponential in propagation distance and frequency dependent. When the 
generalized model is applied to ocean wave statistics, the data are 
scattered within such a small band of the frequency scale that func- 
tions supporting almost all proposed theories can be fitted equally well 
to the data. 

Accordingly, the primary objective of this research was the collec- 
tion of attenuation data in a small scale real environment where the 
measurable frequencies are considerably displaced from those normally 
measured in ocean swell in order to give better definition to the at- 
tenuation function. This requirement for collection of data gave rise 
to the second objective, the development and evaluation of a system for 
recording and analyzing waves in a small scale real environment. 
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II* CONCEPTS 



A. COMPARISON OF ENVIRONMENTS 

Ocean wave attenuation studies involve the measurement of the energy 
of the relatively low frequency components of swell propagated over long 
distances from the area of generation. It is known that these low 
frequencies required high winds of appreciable duration for their develop- 
ment, hence the source of the swell waves can be located with fair ac- 
curacy in relation to the distance from the measuring instrument. A 
knowledge of wind or pressure conditions in the fetch can provide a 
reasonable estimate of the energy of waves in the fetch. 

High frequency waves can be generated by the prevailing low veloc- 
ity winds which are normally found over most of the ocean surface. For 
this and other reasons of distance and instrumentation, it would be ex- 
tremely difficult to make meaningful measurements of these frequencies 
with sufficient precision in the ocean to determine attenuation. For 
purposes of this paper, high frequencies shall be defined as those fre- 
quencies approaching, but not including, frequencies in which surface 
tension is the predominant restoring force rather than gravity. 

A lake or non-tidal estuary of moderate depth provides an ideal 
facility where experimental control is possible in the comparative 
measurement of high frequency deep water waves, although the scale of 
wind circulation in relation to the size of the lake or estuary provides 
no geographical separation between areas of generation and non-genera- 
tive propagation. On this scale, however, the mean wind may change 
rapidly in speed and/or direction. In a case where the wind speed 
diminishes rapidly, some frequencies generated at the higher wind speed 
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would receive no more energy and would propagate as swell. This con- 
cept, that fetch and swell can be separated by time just as there is 
the geographical separation in the ocean, is a fundamental assumption 
of this paper. 

B. THE GENERALIZED ATTENUATION FUNCTION 

Let S(f,x) represent the energy density spectrum computed from a 
wave record at a downwind (propagation) position, and let S(f,0) represent 
the energy density spectrum computed from a wave record at an upwind 
(generation) position, where f is the component frequency and x is the 
separation between two points of measurement along the axis of maximum 
propagation (the direction of the mean wind during generation). Then, 

S(f,x) « S(f ,0) k(x,0 , f ) G(f,x) (1) 



where G(f,x) is the attenuation function, and k(x,0,f) denotes what 
portion of S(f,0) would reach the downwind position if no attenuation 
were involved. This is a function of geometry and the beam pattern in 
the generating area. In this paper the model of Pierson, Neumann and 
James (9) will be used, where for all f 



k(x,0,f) 



_ 2 _ 

TT 




0 d0 



( 2 ) 



The angles 0^ and 0^ are measured at the edges of the fetch from the 
direction of the axis described above to the direction of the downwind 
measurement point. 

Let P(f,x) be a prediction of the spectrum at the downwind position 
which has not allowed for attenuation. Then, 



P(f,x) = S( f ,0) k(x,e,f) H(f,x) 



(3) 
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where H(f,x) is an error function accounting for possible statistical 
errors in measurement at the upwind and downwind positions, and in- 
cludes differences between the angular spreading model and the actual 
angular spreading phenomenon. The model has proved to be fairly representa- 
tive of actual geometrical spreading, and H(f,x) should be at most a very 
slowly changing function of x. The spacial variation of H is thus ignored. 

Divide (3) by (1) 



R(f *x) 



P(f.x) 

S(f ,x) 



H(f ,x) 
G(f ,x) 



(4) 



Differentiate 



dR 

dx 



-H(f,x) 



G'(f,x) 

G 2 (f,x) 



Substitute (4) in (5) 

= -R(£,x) 



Now, assume that R(f,x) can be represented by: 

r»/r \ r. mX 

R(f,x) = C e 



(5) 



( 6 ) 



(7) 



where m - m(f), and C is a constant. This assumption appears to be 
valid from empirical data in ocean swell studies over a wide range of x. 
Differentiating (7) and substituting (7) in the result, 

R*(f,x) - Cm e mX = m R(f,x) (8) 

then substitute (8) in (6), and 

G * (f ,x) = -m G(f,x) (9) 



17 



which is a simple differential equation with solution 



G(f,x) = G Q (f,0) e" mX (10) 

Clearly, at the upwind measuring position, the function G o (f,0) = 1 
and 

G(f,x) = e' mX (11) 

The attenuation coefficient m is a function of frequency and environ- 
mental conditions. Empirical evidence indicates that it increases with 
frequency, and from physical considerations it should reasonably in- 
crease from following wind to opposing wind conditions. Other environ- 
mental conditions may also be important. 



18 



Ill, INSTRUMENTATION 



Two resistance wave meters 9 adapted from a design by Pearlman (5), 
were constructed and used to collect data for this paper. A schematic 
diagram of the circuit used in both meters is shown in Figure 1. 

Pearlman J s paper demonstrated that this type of wave sensor is es- 
sentially linear over a range of wave heights of a few inches in static 
and vertical oscillation tests and has errors of less than one percent 
when tested in orbital motion. Static tests of the modified design used 
in this research confirmed the linear characteristics of the instruments 

The fixed resistors in the wheatstone bridge were standard 57 0 toler 
ance resistors* but were matched by measurement on a resistance bridge 
to have values within 0.5 ohm of the resistance shown. The variable 
resistor was used to null the bridge at an initial immersion in the 
field prior to recording data. 

The signal generators used were Hewlett Packard Model 2000 audio 
oscillators. The recorder was a 14 channel Ampex CP-100 tape recorder. 
As the signal generators and the tape recorder both have one grounded 
terminal, the isolation transformer was required to prevent grounding 
one leg of the wheatstone bridge. 

Alternating current was used to minimize electrolytic effects on 
the electrodes. The full wave rectifier shown in Figure 1 demodulated 
this carrier frequency, passing a direct current signal to the re- 
corder that varied in response to the waves. The signal generators 
were operated at 3 volts (RMS) at about 2000 Hz; however, when operat- 
ing two wave meters on one recorder, it was necessary to have about 150 
Hz difference in the output of the signal generators to prevent a beat 
frequency from interfering with the wave records. 
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FIGURE I. SCHEMATIC OF WAVE METER CIRCUIT 



A dual beam oscilloscope was used to monitor the input and output 
of the tape recorder to insure proper recording of the wave signal. The 1 
Mfd capacitor used across the rectifier output matched the input impedance 
of the tape recorder and oscilloscope, but the use of other equipments 
with different input impedances would require a change in capacitance to 
provide a smooth rectifier output without attenuating wave frequencies. 

The Ampex CP-100 tape recorder was operated with frequency modula- 
tion recording and reproducing amplifiers for the wave record channels 
at a tape drive speed of 1-7/8 inches per second. A third channel was 
used with amplitude modulation (direct) amplifiers for a voice time 
record. 

Details of the sensing probes are shown in Figure 2. The frames 
were designed to provide adequate rigidity, yet to minimize turbulence 
and reflections. The lower horizontal supports were designed to be at 
a depth greater than one-half wavelength of the anticipated waves. The 
fishing line shown was used as an insulator and to minimize the amount 
of electrode wire immersed. If the electrodes had extended to the lower 
horizontal support, the wheatstone bridge could not have been brought 
to near-balance by the variable resistor. In practice, it was found 
that the bridge could be nulled at an electrode immersion of no more 
than 2 centimeters. 

The probe frames were attached to 1 inch diameter galvanized pipes 
driven into the lake bottom by means of the 1 \ inch diameter sleeves 
shown in Figure 2. The meter circuits were mounted on top of the frame 
and connected to the oscillators and recorder ashore with Alpha #2311 
low capacitance shielded microphone cable laid on the bottom. 



21 




22 



IV. DATA COLLECTION 



Roberts Lake, 1500 feet from Monterey Bay in Seaside, California, 
was used for the collection of data in support of this paper. The lake, 
shown in Figure 3, is uniformly 5 to 5% feet deep, with a soft mud bottom. 
Dense reed growth in the water around the circumference minimizes wave 
reflections. 

The Monterey Bay area was considered suitable for this small scale 
attenuation study because its climatological characteristics are favor- 
able to sudden changes in the wind. In addition to unpredictable changes 
in the wind field, diurnal heating and cooling of the Salinas Valley 20 
miles inland results in a sea breeze during the day which dies suddenly 
at dusk. This general characteristic is modified or not apparent when 
the central California coast is under the influence of transient high or 
low pressure regions. In three weeks of daily attempts to collect data, 
only one day produced the desired rapidly dying wind. 

Two sets of data, identified as data sets A and B, were collected 
on 21 November 1968. The wind speeds, as measured with a Casella pre- 
cision anemometer mounted eight feet above the lake surface, are shown 
in Figures 4 and 5. The Casella anemometer provides a mean wind speed 
by recording the revolutions of a three cup rotor on a digital counter. 

The wind data in data set A, Figure 4, were recorded for one-minute aver- 
ages, and in data set B, Figure 5, for half-minute averages. 

Details of wave sensor placement are shown in Figure 3. For data 
set A, the sensing probes were on an axis of 230° (M) . Although the 
actual distance between probes was 400 feet, the direction of the wind 
(from 300°(M)) reduced the effective propagation distance to 138 feet. 

For data set B, the probes were 120 feet apart on an axis of 290° (M). 
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Prior to collecting the wave data, both probes were immersed 
simultaneously in relatively still water near the edge of the lake at 
two levels of immersion and one minute recordings of the signals were 
placed on magnetic tape for calibration of the wave meters. The average 
values of the one minute calibration signals were later used to scale 
the wave records to actual amplitudes in millimeters. 
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FIGURE 4. WIND RECORD FOR DATA SET 
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FIGURE 5. WIND RECORD FOR DATA SET 



V. ANALYSIS OF DATA 



A. DIGITIZING 

The calibration records and selected portions of the wave records 
were digitized with an SDS 9300-COMCQR 5000 hybrid computer system. 

The Ampex CP- 100 tape recorder was patched to the C0MC0R 5000 analog 
computer through operational amplifiers included in the logic board of 
the computer and a low pass filter constructed from logic board components. 
The amplifiers were adjusted to amplify the maximum signal on the magnetic 
tape record to slightly less than the maximum input voltage of the analog 
computer. 

A schematic diagram and plot of the measured characteristics of the 
low pass filter are shown in Figure 6, To minimize aliasing in later 
spectral calculations, 128 Hz was selected as a filter cut-off frequency. 
Although frequencies of that order were far beyond the range of expected 
surface wave frequencies, noise could possibly have existed in that range 
from either inherent equipment noise or beats between the wave meter 
oscillators. 

To provide a minimum of two samples per cycle up to the cut-off 
frequency, the record was sampled at a rate of 256 samples per second. 

The digital computer sampling program used was developed by Post (10). 

B. CONVERSION OF DATA FOR IBM 360 COMPUTER SYSTEM COMPATIBILITY 

The core capacity of the SDS 9300 computer used for digitizing the 
data was far too small for spectral analysis of the wave records (see 
next section). Therefore, use was made of the much larger core capacity 
of the available IBM 360 computer system. The two computers do not share 
common symbology and storage parameters; however, a program developed by 
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FIGURE 6. ANALOG FILTER CHARACTERISTIC 
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Post (10) for conversion of data from one system to the other was used 
to transfer and convert the digitized data from the tape record made on 
the SDS 9300 computer to a tape record compatible with the IBM 360 
computer system. 



C. CALCULATION OF ENERGY DENSITY SPECTRA 

In a particular finite segment of a wave record of length let 

A(t^) represent the time series of N discrete wave amplitudes which are 
spaced at intervals of time At. Then the time of each sample amplitude 
is 

t^ = k At, k = 1 , . . . ,N (12) 



The record is transformed from the time domain to the frequency 
domain by the Fourier transform 



^( f j> " TT 



N-l 

I 

k = 0 



A(t k ) 



-2 TT i jk/N 



,J-0,...,N-1 (13) 



where 



and 



f . = j Af are the component frequencies 

^(f.) are the complex amplitudes of the f. 
J J 



In the computer, this transform is accomplished by the IBM library 

subroutine HARM, which uses the fast Fourier transform algorithm developed 

by Cooley and Tukey (3) 9 and requires that N be an integral power of 2, 

^(f^) is a two-sided amplitude spectrum where N/2 corresponds to the 

aliasing or Nyquist folding frequency, f , which is the highest frequency 

c 

analyzed. If the wave record had component frequencies above f ^ , the ampli- 
tudes of these frequencies would be aliased into the spectrum ft(f^). A 
Nyquist frequency of 128 Hz was selected on the basis of the analog filter 
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characteristics as discussed in section A. The Nyquist frequency deter- 



mines the minimum sampling interval, At, by the relation 

At = l/2f 

c 

The energy density spectrum, S(f is given by 

S(f j> ■ 2 !*»!> , J - 0 N/2 

A f 



(14) 



(15) 



where is the complex conjugate of ^(f j) ar *d the factor 2 accounts 

for the inclusion of only half of the two-sided symmetrical amplitude 
spectrum. The S(f .) represent the average energy per unit frequency in 
an elementary frequency band of width Af centered on f.. 

The S(f^) one realization of an ergodic random process is a set 
of statistics, each of which is considered to have the Chi-square distri- 
bution with two degrees of freedom. The number of degrees of freedom 
determines the confidence limits for a statistic, and a meaningful compari- 
son of statistics requires that the range of confidence limits be suffi- 
ciently small. The relation between the number of degrees of freedom 
and the confidence limits for the Chi-square distribution is given by 
Pierson and Marks (8). 

An accepted method of increasing the confidence in the statistical 
estimate of energy density is to average a number, M, of elementary 
frequency bands into a wider equivalent band of width B^, When this is 
done, each bandwidth 8^ has 2M degrees of freedom. As the two degrees 
of freedom in each elementary bandwidth reflect two independent samples 
in the time domain, the record length, T , required for the analysis is 
then specified by: 



T = NAt = 2 f M A t /B 

r c e 



(16) 
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Obviously, the increased confidence is bought at the expense of defini- 
tion,, 

The parameters used for analysis of the wave records collected on 
21 November 1968 are summarized as follows: 

T =64 seconds 
r 

N = 16384 

£t = 1/256 seconds 

f = 128 Hz 
c 

B = 0.5 Hz 
e 

M « 32 

For 64 degrees of freedom, the 957> confidence limits are 1.5 and 
0.7, which means that for an infinite number of realizations of this 
particular wave process, 95? 0 of the realizations of the spectrum would 

fall between 1.5 S(f.) and 0.7 S{f.). 

J J 

In the discussion of the results of the spectral analyses in the 
following section, confidence limits have been omitted for clarity, but 
their existence must be recognized. In fact, the computed S(f^) are t * ie 
most likely values. 
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VI. RESULTS AND IMPLICATIONS OF SPECTRAL CALCULATIONS 



In the text and figures of the discussion that follows, the time 
(hours :minutes : seconds) assigned to a spectrum is the mid- time of the 
64-second wave record from which the spectrum was derived. For the wind 
data subset of data sets A and B, refer to Figures 4 and 5, respectively. 

A. COMPARISON OF SIMULTANEOUS UPWIND AND DOWNWIND SPECTRA 

Figure 7 shows the energy density spectra computed from simultan- 
eous realizations at the upwind and downwind wave probe positions at 
11:55:56 and 11:57:00* This comparison implies that spectra from down- 
wind probe data may be underestimated by a factor of approximately two. 
The most obvious cause for this under-estimation is an error in calibra- 
tion* however, it is possible that the difference may be an effect of 
the geometry of the lake. 

B. ATTENUATION INDICATED BY DATA SET A 

Figure 8 shows the energy density spectra at the upwind probe posi- 
tion at 11:57:00 when the average wind at the eight foot level had been 
four knots or more for at least 30 minutes, and at the downwind probe 
position 112, 176, and 240 seconds later. During the four-minute inter- 
val represented, the wind diminished to 1.3 knots. 

In the development of the generalized attenuation function earlier 
in this paper, the true energy density spectra at the generation and 
propagation positions were related by: 

S(f,x) = S(f,0) k(x,0,f ) G(f,x) (1) 

and the attenuation function was specified by: 

G(f,x) = e' mX (11) 
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FIGURE 7. COMPARISON OF SIMULTANEOUS SPECTRA 
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ENERGY DENSITY, MM Z -SEC 




FIGURE 8. ENERGY DENSITY SPECTRA- DATA SET A 
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If it: is assumed that the energy density spectra, S(f .,0) and 
S(fj,x), derived from measurements of the upwind and downwind probe 
positions, respectively, are equivalent to the continuous spectra, then 



G(f »x) = 



s(r,x) 



S(fj»0) k(x s ,9 9 £) 



= e 



-mx 



(17) 



which, as the attenuation coefficient, m ~ m(f), is sought, can be 
written s 



In 



m(f j ) = 



S(f .,0) k(x ,9 , f ) 

~ S(fj,x) 



(IB) 



A valid comparison of energy densities for the determination of 
m(fj) can only be made when the identical wave components are present 
at both measurement points. Therefore, S(f^,0) measured at the upwind 
probe can only be compared with S(fj,x) measured at the downwind probe 
after the time, t, required for the energy of the component frequency 
f to propagate over the distance x, the effective distance between probes. 



This is given for gravity waves by 



t = x / V = 4 x 7 / f / g (19) 

S 

where V is the group velocity of the f. and g is the acceleration of 
g J 

gravity. 

Table 1 summarizes the data and computations for m(f ^) f° r data set 
A. As the S(f.) were computed at specific times, the frequencies whose 
propagation times corresponded to the time intervals between upwind and 
downwind spectra were computed by the inverse of equation (19) and 
rounded to the nearest % Hz to correspond to the center frequency of the 
band represented by f „ , unless the calculated frequency fell midway be- 
tween center frequencies, then rounding was done in both directions and 
dual calculations made . 
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The second calculation of m(f.) under each frequency uses 2 S(f.,x) in 




wind spectra mentioned in section A, 



TABLE 1 



SUMMARY OF STATISTICS AND ATTENUATION COEFFICIENTS FOR DATA SET A 







t 


X 


f . 
J 


S(f.,0) 




S(f j,x) 


m( f j) 


Upwind 

time 


Downwind 

time 


sec 


ft 


-1 

sec 


mm^/Hz 


k(x,0. 


f ) 2 
mm /Hz 


ft' 1 


11:57:00 


11:58:52 


112 


138 


2.0 


1.31 


.95 


.648 


.00473 


















-.000348 


11:57:00 


11:59:56 


176 


138 


3.0 


2.24 


.95 


.136 


.0199 


















.0149 










3.5 


2.35 


.95 


.0349 


.0301 


















.0251 


11:57:00 


12:01:00 


240 


138 


4.5 


.519 


.95 


.00265 


.0377 


















.0327 



Now, as m = m(f), assume that m(f) can be expressed by the equation 

ra(f) = K f P (20) 

where K is an attenuation constant and p is an integer. 

Table 2 shows the values of K for the values of m(f) given in 
Table 1 for p = 2,3,4. 
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TABLE 2 



f 

-1 

sec 


m(f ) 
ft" 1 


ATTENUATION 


CONSTANTS FOR DATA 


SET A 


P = 4 
K 

sec^/nm 


m(f ) 

-1 

nm 


p - 2 

K 

sec^/nm 


p - 3 
K 

sec^/nm 


2.0 


.00473 


28.4 


7.1 


3.55 


1.77 




-.000348 


-2.08 


-.52 


-.26 


-.13 


3.0 


.0199 


119. 


13.4 


4.41 


1.47 




.0149 


89.3 


9.9 


3.31 


1.10 


3.5 


.0301 


180. 


14.7 


4.18 


1.20 




.0251 


150. 


12.2 


3.49 


1.00 


4.5 


.0377 


226 


11.1 


2.46 


.545 




.0327 


196 


9.6 


2.13 


.472 



Co COMPARISON WITH OCEAN SWELL DATA 

Some energy density spectra of recorded Pacific Ocean swell from 
fetches at known distances from the swell observation point were made 
available by Wickham (12) o Values of m and K for these data have been 
computed in the manner prescribed for these same values in Tables 1 and 
2 above, and are displayed in Table 3. The data for Table 3 may include 
situations in which the winds in the propagation areas were across or 
opposed to the direction of propagation, where the data for Tables 1 
and 2 were for only the following wind situation. 
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TABLE 3 



ATTENUATION COEFFICIENTS AND CONSTANTS FOR OCEAN WAVES 



f 


m(f) 


P = 2 


P = 3 


II 






K 


K 


K 


-1 

sec 


-1 

run 


sec^/nm 


sec^/nm 


sec^/nm 


.105 


.00275 


.275 


2.75 


27.5 


.095 


.00224 


.248 


2.62 


27.5 


.095 


.00261 


.289 


3.04 


32.0 


.085 


.00266 


.369 


4.34 


51.0 


.075 


.00378 


.670 


8.95 


120. 


.055 


.00115 


.380 


6.93 


125. 



D. RESULTS OF DATA SET A AND TABLE 3 

1. The possible error of a factor of 2 between spectra as dis- 
cussed in section A. does not have a significant effect on the deter- 
mination of attenuation for data set A. With the exception of the values 
at fj = 2.0 Hz, the attenuation coefficients, m(f) are all of the same 
order of magnitude; and with the exception of the negative value for the 
second estimate of K for 2.0 Hz, the attenuation constants for each power 
of frequency are of the same order of magnitude. 

2. If multiplying S(fj, x) by 2 had a small effect, then variations 
within the confidence limits for the spectral estimates will also have 
small effect. 

3. It is quite evident that if an attenuation function can be 

described over a wide range of frequencies by a simple relationship such 

3 

as proposed by equation (20), then it must be a function of f . 



39 



E. DATA SET B 

The wind in data set B died more slowly than the dramatic change 
which allowed the extended calculations for data set A. For those short 
time periods where the wind declined rapidly, S(fj) were computed and 
are displayed in Figures 9, 10, and lie As each figure represents only 
one attenuation observation, the corresponding values of m(fj) anc * K are 
shown on the figures. 

Less confidence is placed in the results of data set B, as the 
results are more erratic and the spectra seem to be influenced near the 
frequencies of interest by a low frequency disturbance that is not evi- 
dent in data set A. This low frequency could have been due to a surge 
set up in the lake by the gusting of the wind during its decline. 

F. ADDITIONAL OBSERVATIONS 

The persistence of the energy density at 2.0 Hz is interesting, 
especially in view of the negative values obtained when 2 S(fj,x) is 
substituted for S(fj,x) for the calculation of m(fj) and K. These 
negative values indicate that the energy in that particular frequency 
may be increasing with distance. A possible cause for this is that 2 Hz 
may be a low enough frequency to suffer negligible attenuation in such 
a short distance. 

The second peak that occurs in both data sets at 6 Hz is also note- 
worthy. Possible explanations for this are that a frequency of 6 Hz is 
beginning to be influenced by surface tension effects, or the effect may 
be a result of non-linear interactions as have been reported by Pierson 

(7). 
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FIGURE 9. ENERGY DENSITY SPECTRA - DATA SET B 



u 
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FIGURE 10. ENERGY DENSITY SPECTRA -DATA SET B 



42 



ENERGY DENSITY. MM - SEC 
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FIGURE II. ENERGY DENSITY SPECTRA— DATA SET B 
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The decision to establish the Nyquist frequency at 128 Hz was not 
unwarranted, as examination of the spectra beyond the water surface wave 
frequencies show occasions of significant energy, particularly around 
60 Hz. 



44 



VII* CONCLUSIONS AND RECOMMENDATIONS 



The instrumentation, recording, filtering, digitizing and analyti- 
cal techniques employed in support of this investigation were extremely 
effective in providing well defined energy density spectra from waves in 
a small scale environment. 

The dying wind situation gives a definable separation of wind wave 
generation and swell propagation for high frequency waves when the wind 
dies rapidly and smoothly. 

Strong evidence is provided in support of a contention that an 
empirical attenuation function can be applied to a wide range of fre- 
quencies, and it is proposed that this function might be 



G(f,x) = e 



-Kf 3 x 



where K is approximately 3 sec /nm for following winds and may increase 

3 

to as much as 8 sec /nm for opposing winds. Future extensive investi- 
gation may show the frequency exponent to not be an integer, but for the 
current state of definition of the attenuation function this integral 
representation is satisfactory. 

It is recommended that further studies be undertaken to provide 
better definition of such an empirical attenuation function, with emphasis 
placed on intermediate frequencies between the high frequencies investi- 
gated in this paper and the lower frequencies more frequently studied in 
the ocean. 
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